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Nonequilibrium dynamics of a Lieb-Liniger system in the presence of the hard-wall potential is 
studied. We demonstrate that a time-dependent wave function, which describes quantum dynamics 
of a Lieb-Liniger wave packet comprised of N particles, can be found by solving an iV-dimensional 
Fourier transform; this follows from the symmetry properties of the many-body eigenstates in the 
presence of the hard-wall potential. The presented formalism is employed to numerically calculate 
reflection of a few-body wave packet from the hard wall for various interaction strengths and incident 
momenta. 
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I. INTRODUCTION 

Recent intensive theoretical studies of one-dimensional (ID) Bose gases have been motivated mainly by the experi- 
mental realizations [H, 0, H, B 13 of the models [H, 0] describing these systems. Ultracold atomic gases in tight atomic 
wave guides can have their transversal degrees of freedom essentially frozen such that their motion becomes effectively 
one-dimensional P, S S, H Q. These systems are described in terms of the Lieb-Liniger model, i.e., a system of 
identical Bose particles in ID with pointlike ^-function interactions of arbitrary strength c Q. In the limit of strong 
coupling (c — ► oo), the Lieb-Liniger gas approaches the Tonks- Girardeau regime of impenetrable bosons 0], which 
has been experimentally achieved as well |2j; the Tonks- Girardeau regime occurs at very low temperatures, low linear 
densities, and with strong effective interactions [illsjicj. Experiments are even capable of exploring nonequilibrium 
quantum dynamics of these ID many-body systems 3, |4[, which may occur after some sudden change in the system's 
parameters. These ultracold atomic assemblies are well isolated from the environment, that is, their quantum coher- 
ence stays preserved for long times. Therefore, they may serve as a playground to investigate relaxation of isolated 
quantum many-body systems, which is one of the most interesting questions in theoretical physics (e.g. see Refs. 
jlH H2I HH , [lj, [TH, [TH, [I?], EH an d references therein). Subsequent relaxation of ID gases via collisions is greatly 
determined by the reduced dimensionality and the integrability of the underlying models. We are motivated to study 
the time-dependent Lieb-Liniger model because (i) today's experiments can explore fundamental physical questions in 
these systems [H, 0| , and (ii) one can construct exact solutions of some relevant problems for all interaction strengths 
(from the mean field regime up to the strongly correlated regime) [ljl [2(| [22|, [H, [24| . 

The eigenstates of the Lieb-Liniger model (without an external potential present), which were constructed by 
employing the Bethe ansatz Q , are determined by a set of quasimomenta; when periodic Q boundary conditions arc 
imposed, the quasimomenta must obey a set of transcendental Bethe equations @ . The Lieb-Liniger eigenstates in the 
presence of the hard- wall (i.e., on the semi-infinite line) can be constructed via superposition of free space eigenstates 



[251 ] ; again, if the quasimomenta should obey a particular set of transcendental equations [25j , this superposition yields 
eigenstates in an infinitely deep box 12511 . Recent years have witnessed an increasing interest in exact solutions of 
these models (e.g., see (26l l27l. l28l. [29l [30L l3lL l32j and references therein), most of which are focused on the properties 
of the ground and excited eigenstates (see also Refs. [33|,|34|). Unfortunately, the Lieb-Liniger model does not reveal 
exact solutions in the presence of some external trapping potential V(x) (e.g., the harmonic potential). 

In the Tonks- Girardeau limit c —* 00, the methods for finding eigenstates Q, time-dependent solutions [35|, as well 
as observables (e.g., see Ref. [3(| for the system of hard-core bosons on the lattice and |37j for the continuous Tonks- 
Girardeau model) are much simpler due to the Fermi-Bose mapping, which in a simple fashion maps a fermionic wave 
function describing spinless noninteracting fermions onto a Tonks- Girardeau wave function fjj, |35|. It is important 
to emphasize that these methods are valid for any external potential. Perhaps the simplicity of the methods and 
phenomenological relevance of the model [|| have lead to increasing interest in quantum many-body dynamics of 
Tonks- Girardeau gases. Some of these studies include dynamics during free expansion [36l. l38l. l39l. |40|| . dynamics of 
dark soliton-like states [3l| , and reflections from a periodic potential [37| . 
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In the case of finite interaction strength c, it is far more difficult to calculate exact many-body wave functions 
and/or observables describing dynamics of time-dependent Lieb-Liniger wave packets. Without attempting to pro- 
vide a review, let us mention a few approaches utilized to study noncqulibrium dynamics of ID interacting Bosc 
gases. The hydrodynamic formalism [1 01 ] (the local density approximation) can be formulated in terms of the Non- 
linear Schrodinger like equation with variable nonlinearity [4l| : this approach reduces to the Gross-Pitaevskii theory 
in the weakly interacting limit [HI EH- More sophisticated numerical approaches include the time-evolving block 
decimation algorithm [42[, which has recently been utilized to study relaxation following a quench in a ID Bose 
gas [43| , the twoparticle irreducible (2PI) effective action approach 0, HH , the multiconfigurational time-dependent 
Hartree method for bosons (MCTDHB) [46[ (the MCTDHB method is numerically exact when sufficiently many 
time-dependent orbitals are taken into account), and the multiconfigurational time-dependent Hartree method (e.g., 
see Ref. [I?) and references therein). Reference [24[ provides a discussion of several methods which can be used to 
describe nonequilibrium dynamics of Lieb-Liniger gases with greater focus on the form- factor approach [2~i| , which has 
been recently utilized to calculate equilibrium correlation functions of a ID Bose gas (see (48| and references therein). 
A broader review discussing many-body physics with ultracold gases can be found in Ref. (4j|. We also mention a 
recent review on quantum transients [501 ] . 

An interesting exact method has been outlined by Gaudin way back in 1983 [l9| : A time-dependent Lieb-Liniger 
wave function on an infinite line, in the absence of an external potential, can be constructed by acting with a 
differential operator (which contains the interaction str eng th parameter c) onto a time-dependent wave function 
describing noninteracting (spin polarized) ID fermions [l9l l2ll. I22I [23l . [5l| . For dynamics of a Lieb-Liniger wave 
packet comprised of N particles, this method reduces to finding an A-dimensional Fourier transform, which can be 
used to extract the asymptotic behavior of the wave function and some observables during the course of ID free 
expansion [22|, [23[ . In this article we investigate the possibility of extending this approach to study dynamics of a 
Lieb-Liniger wave packet in the presence of the hard-wall potential. 

Our interest in quantum dynamics in the presence of the hard- wall potential is in part motivated by experiments. 
More specifically, the interaction of Bose-Einstein condensates (BEC) with surfaces is of interest for implementations 
of atom interferometry on chips (52j. A BEC falling under gravity, and then reflecting from a light-sheet, has 
been experimentally and theoretically studied in Ref. [53l ]. Moreover, one of the prominent experimental activities 
nowadays is deceleration of atomic beams by reflection from a moving mirror. This work first started with neutrons 
being cooled by reflecting from a moving Ni surface [54] ]. In cold atoms physics, there have been several experiments 
for manipulation and slowing down atomic beams with the use of reflection mirrors [55l . IBlil . H3] • 

Here we explore, by using exact methods, dynamics of Lieb-Liniger wave packets in the presence of the hard-wall 
potential, more specifically, reflection of a Lieb-Liniger wave packet from such a wall. The outline of the paper is as 
follows. In Sec. |TT]we introduce the model and outline the construction of eigenstates in the given external potential. 
In Sec. IIIII we analytically discuss time-dependent quantum dynamics of the system which starts from a general initial 
condition. By employing the symmetries of the Lieb-Liniger eigenstates, we demonstrate that a time-dependent 
Lieb-Liniger wave packet reflecting from the wall can be calculated by solving an A-dimensional Fourier transform, 
where N is the number of particles. This opens the way to calculate the asymptotic properties of the wave packet 
by employing the stationary phase approximation as in Refs. [22I |23| for free expansion. In Sec. IIVI we utilize the 
formalism to numerically study dynamics of single-particle density and momentum distribution of a few-body wave 
packet reflecting from the wall. We find that the wave packets for smaller interaction strength c get reflected at a 
slower rate, because they get compressed more strongly as the wave packet hits the wall. The interference fringes 
which occur during the dynamics have larger visibility for smaller values of c. 



II. EIGENSTATES IN THE PRESENCE OF THE HARD- WALL POTENTIAL 



The Lieb-Liniger model describes N identical bosons in one spatial dimension, which interact via a repulsive 5- 
function potential of strength c > 0. The model can be represented in terms of the Schrodinger equation: 

i— =-^2q- 2+ X] 2c $( x i - x j) 1 P + V(xi)ip(xi, X N , t). (1) 

i=l i l<i<j<N i=l 

As we have already stated, this model can be experimentally realized with ultracold atoms in tight atomic waveguides. 
The spatial and temporal coordinates (x and t, respectively), as well as the potential V(x) are dimensionless in this 
paper. Their connection to physical units is as follows: x = X/Xq, t = T/Tq, and V(x) = U(X)/Eo, where X, T and 
U(X) are space, time, and energy variables in physical units. Given the mass of the atoms m, the choice of an arbitrary 
length scale X sets the time scale T = 2mJfg//i, and energy scale E — fi 2 /(2mX 2 ). Suppose that the transverse 
confinement of the atomic waveguide is described by a harmonic oscillator with frequency loj_. The interaction 
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parameter c is proportional to the effective ID coupling strength gm [U, 2c = gm/ (XoEq) = gm 2mXo/h 2 , which is 
related to the ID scattering length am via gm — —2Y?/mam\ the ID scattering length am = — {l\/a)(l — Ca/V2l±) 
depends on three-dimensional scattering length a and the transverse oscillator width l± = ^/h/mui± (the constant 
C= 1.4603...). 

In the present paper, we focus ourselves on the dynamics (in time) of a Lieb-Liniger wave packet in the presence of 
the hard-wall potential, that is, 



V(x) 



0, if x > 
co, if x < 0. 



(2) 



We will show that the solution of this problem can be constructed by solving an .ZV-dimensional Fourier transform. To 
this end, we need eigenstates of a Lieb-Liniger gas in the hard- wall potential. First, let us write down the Lieb-Liniger 
cigenstates in free space (i.e., x G (—00,00) without external potentials and any boundary conditions): 



v^> = AA({fc»5>i) p n 



8gn(xj - Xi) + - {k P j - kpi) 



J. H j kpjXj 



= AA({fc})£(-l)MP,{fc})e^ 



(3) 



where {k} = {k rn | m = 1, . . . , ./V} is a set of (real) distinct quasimomcnta which uniquely determine the cigenstate, P 
denotes a permutation of N numbers, P G Sn, and we have implicitly defined a(P, {k}). The normalization of these 
eigenstates is given by [H, [H[ 
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M(2tt) n Y[ 









Af({k}) \ 

that is, within the fundamental sector in fc-space, k\ < • • • < fcjv and k[ < • • • < k' N , we have 



/ ^{k}^{k'}dxi ■ ■ ■ dx N = Y[ 3(kj - k'j 
J-00 J=1 



(4) 



The Lieb-Liniger eigenstates in the presence of the hard- wall (denoted by 4>{k}) were first constructed by Gaudin 
[25l ] as a superposition of 2 N free-space eigenstates. This superposition obeys the hard-wall boundary condition: 
</ ) {fc}(2 ; i = 0, £2, . . • ,xn) = in the fundamental sector Ri : x\ < x% < . . . < xn of a;-space. These eigenstates are 
expressed as follows: 



^{fc} =]L A ({ £ }'W)^}' 



(5) 



where {e} = {e m \ e m G { — 1, 1}, m = 1, . . . , N} and {ek} = {e rn k m \ e m G {—1, 1}, m = 1, . . . , N}] evidently, there are 
2 N such sets and therefore 2 N terms in the sum ([5]). The quantity A({e}, {k}) is defined by 



where 



A({e}, {k}) = ei • ■ • e N A'(e 1 ki,e 2 k 2 , £jvfcjv)i 



A'(k l7 k 2} . . . ,k N ) 



in 



kj+ki 



(6) 



(7) 



are the coefficients utilized in the superposition. It is straightforward to verify that indeed 4>{k}( x i = 0) x 2, ■ ■ ■ , xn) = 
in the fundamental sector i?i [2{|. 

However, it is not simple to prove that these eigenstates are orthogonal and normalized. This is of key importance 
if one wishes to project some initial state onto these eigenstates and calculate time-evolution in the standard fashion 
via superposition over eigenstates. In Section llVI we discuss the normalization of eigenstates ((5|), and based on our 
numerical investigations conjecture that these eigenstates are orthogonal and normalized. 
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III. MANY-BODY DYNAMICS IN TIME VIA A FOURIER TRANSFORM 



In this section we demonstrate that a solution of the time-dependent equation {T]) with the hard-wall potential ([2]) 
can be expressed in terms of an iV-dimensional Fourier transform. We assume that at time t = the wave packet is 
localized in the vicinity of the wall. For example, the initial state tpQ can be the ground state wave function in some 
external trapping potential; if at t — this potential is suddenly turned off, the wave packet will start expanding and 
some of its components will be reflected from the wall which will give rise to interference effects. Such a scenario is 
possible to create with today's experimental capabilities [3]. One possible (similar) scenario is as follows: suppose 
that at t = the aforementioned trapping potential is turned off, and that in the next instance the many body wave 
packet is given some momentum kick, say towards the wall; the reflection and interference phenomena will depend 
on the interactions and imparted momentum. During the reflection, particles will collide and one may ask to which 
extent will the initial conditions be forgotten (or blurred) after the reflection? 

To describe quantum dynamics from the initial conditions described above, we write the initial state ipo as a 
superposition over complete set of eigenstates <p{k}'- 



ipo = / b(ki, . . . , k N )cp {k} dk 1 . . . dk N . (8) 

Jo<k 1 <-<k N 

The subsequent derivation is based on the following two relations obeyed by the eigenstates 4>{k} : 

<f>{k} = (-lf<f>{Pk}, (9) 

and 

4>{k\ = -0{fc 1 ,...,fc J _i,-fc J ,fc J+ i...,fc ]V }- (io) 

Equation ^ follows from the definition of (p{ k } in Eq. ([5]), and the fact that the Lieb-Liniger eigenstates in free space 
^{fc} obey identical relation: ip{ k } = ( — l^V'fPfc}! this identity can be traced to the fact that ip{ k y are antisymmetric 
with respect to the interchange of any two variables ki and kj The derivation of Eq. (|10[) is straightforward. Let 
us define a set {e'}, which corresponds to the set {e} as follows: {e' 1; . . . , e^y} = {ei, . . . , ej_i, — ej, ej-^i, . . . , e^r}; 
it is evident from the definition ([5]) that = J2{e>} {^})"0{e'fc} ■ Furthermore, let us denote {k'} = 

{fci, . . . , kj—i, —kj, kj+i . . . , fc_/v}, i.e., the set of fc-values {k'} is identical to the set {k} except that kj is reversed 
in sign. By using A({e}, {k}) = — A({e'}, {k'}) and {e'fc'} = {ek} we have 

hk> } = £ A({e'}, {fc'})w fc ' } = - E A ({ £ >' = -hk}' (ii) 



that is, we obtain Eq. (JTUJ) . We note in passing that if any kj = 0, then <f>{ k } = 0, which follows from Eq. (JTI 
furthermore, <p{k} IS a ls° zero whenever any two of the quasimomenta fcj and kj are equal. 

Due to the symmetry of the hard- wall eigenstates 4>{k} presented in Eqs. © and (jTOJ) , a complete set of eigenstates 
is spanned in the region of the fc-spacc defined by < ki < ... < fcjv , which we will refer to as the fundamental 
region in fc-space, and denote it with . Hence, the integral in Eq. ([5]) spans over Qf . Furthermore, by employing 
relations (O and (jTUJ) , ipo can be written as an integral over the whole fc-space: 

/oo 
dfci • • ■ dk N G{kx, . . . , k N )e l ^ k ^ , (12) 
- OO 

where the function G is defined as 

G(k 1: ...,k N ) = b(h, . . . , k N )A'(kx,. . . , k N )Af({k}) Y[ 



1 

sgn(xj - Xi) + - {k 



(13) 



From Eq. lfl3l) it immediately follows that the time-evolution of a Lieb-Liniger wave packet in the presence of the 
hard-wall can be calculated from an TV-dimensional Fourier transform. 

In order to derive Eqs. (|12jl and (|13[). first note that due to (J9j> and (jlOjl . the projection coefficients satisfy 

b(h, ...,k N ) = {-lfb{k P1 , k PN ), (14) 

and 



b(ki, . . . , k N ) = -b(ki, . . . , kj-i,-kj, kj+i, . . . , k N ); 



(15) 
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the latter identity can conveniently be rewritten as 

b(ki, . . . , k N ) =ei--- e N b{e 1 k 1 , . . . , e N k N ). (16) 

By employing the symmetries of the Lieb-Liniger hard-wall eigenstates, which are inherited by the expansion 
coefficients b(kx, . . . , fejv)j Eq. © can be rewritten as follows: 

^° = Tn I dkl '" dk N b{k x , ...,k N ) V A({e} 7 {k}) 

iV - t /fc 1 >0 1 ...,fc JV >0 | £ j 

x Af{e 1 k 1 , e N k N ) J2 (-l) P a(^ Uk})e 1 ^ (17) 



1 

JV ■ ■'k 1 >0,...,k N >0 



dki--- dk N \ b(exkx, e N k N )A'(exk u e N k N ) 



M 



x Af(exkx, e N k N ) £ (-l) p a(P, {efc})e 4 ^ (18) 



1 

TV! 



x 



dki ■ ■ ■ dkpj b(kx, ■ ■ ■ , kTf)A'(ki, . . • , fc^r) 
. • • , fciv) V (-l) P a(P, {fc})e 4 ^ fep ^ (19) 



p 



1 f 00 

= W\ / *i • • • *JV 6(fcpi, ■ ■ ■ , fcpjv)^4'(fcpi, ■ ■ ■ , fcpjv) 

■'»• p ./-oo 

x AA(fc P1 , . . . , fc PJV )a(P, {fc})e^ fe "^ (20) 
= AH / * P1 ' ' ' dkpN fe ( fc -Pi' ■ ■ ■ ' fcpAf)-4'(fcpi, • • • , k PN ) 

p J -oo 

x AA(fc P1 , . . . , k PN )a(P, {k})e 1 ^ kp * x *, (21) 



from which we immediately obtain Eqs. (fl2|) and ([T3[) because the sum over all permutations P is a sum over AT! 
identical integrals. In the derivation above, the first identity, Eq. (fT7|) . follows from the properties ((9]) and (JT4J). The 
second identity (fT8|) is due to (fPoj) and the definition of A({e}, {fc}) in Eq. ([6]). By employing Eqs. (fT0|) and (Tl6j) . 
the sum over {e} in Eq. fTS]) can be replaced by integrating over the whole fc-space to obtain the third equality, Eq. 
(fT9|) . By using identities A'(kpx, . . . , fcpjv) = A'(kx, ■ ■ ■ , fcjv) and Af(kpx, ■ ■ ■ , ^paO = J\f(kx, ■ ■ ■ , fejy), together with 
Eq. (Till), we obtain ([20]) . 

The time-dependent solution of the many-body Schrodinger Eq. |T]) with V(x) given by @ is simply 

poo 

ij)= / dk 1 ...dk N G(k ll ...,k N )e l ^ {k ' x '-^ t) . (22) 



Thus, by knowing the function G which contains all information about the initial condition, and which is simply 
related to the projection coefficients b(kx, ■ ■ ■ , fciv) of the initial state onto hard-wall Lieb-Liniger eigenstates 4>{k}, 
we can compute the time-dependent Lieb-Liniger wave function in the hard-wall potential by employing the Fourier 
transform. With this identification, an exact analysis of this many-body problem is at least conceptually considerably 
simplified. 

We note that the asymptotic behavior of the many-body state and the observables such as single-particle density 
or momentum distribution can be straightforwardly extracted from expression (|22[) by using the stationary phase 
approximation, as it was done in Refs. [22], l2.'il | for the case of free expansion of a Lieb-Liniger gas [e.g., see Eq. 
(15) in Ref. and Eqs. (18) and (19) in Ref. [H]. From these methods, and Eqs. ([22]) and JT3J), it follows 

that the initial conditions are imprinted into asymptotic states. It is straightforward to infer that the asymptotic 
wave functions, ipooivij • • • i VNi t) = ij)(r)xt, . . . , rjpjt, t) for sufficiently large t, vanish at the hyperplanes of contact 
between particles rji = rjj (i ^ j), which is characteristic for Tonks- Girardeau wave functions [7[. However, it should 
be emphasized that the properties of such asymptotic states can considerably differ from the physical properties of 
a Tonks-Girardeau gas in the ground state of some trapping potential (see also the second item of Ref. [2f|). 
Moreover, the asymptotic momentum distribution coincides, up to a simple scaling transformation, with the shape of 
the asymptotic single-particle density in x-space, reflecting the fact that the dynamics is asymptotically ballistic [23| ; 
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FIG. 1: (color online) Density evolution of a Lieb-Liniger wave packet comprised of N = 3 bosons, which is given some 
momentum kick K (per particle) towards the wall. Insets correspond to the interaction strengths (a) c = 0.25, (b) c = 3, 
and (c) c = 10. The imparted momentum is K = 1. Red dotted lines are for t — 0, black solid lines are for t = 1, and blue 
dashed-lines are for t = 2. 




FIG. 2: (color online) The same as in Fig. [T]but for the momentum distribution 



this means that at asymptotic times, despite of the fact that the wave functions have attained the Tonks-Girardeau 
structure, interactions do not affect the dynamics any more. From the connection between the asymptotic momentum 
distribution and single-particle density one finds that the asymptotic momentum distribution is zero at k = 0, and it 
is located on the positive /c-axis, which simply means that for sufficiently large times the particles move away from 
the wall. 



IV. EXAMPLE: A LIEB-LINIGER WAVE PACKET INCIDENT ON THE HARD WALL 

In this section we study a specific example of a localized Lieb-Liniger wave packet comprised of a N = 3 particles 
reflecting from the hard-wall potential. More specifically, we assume that for t < the Lieb-Liniger system is in the 
ground state of an infinitely deep box denoted by vp g . s .(xi,X2 7 x 3 ). The analytic expression for this ground state was 
found in Ref. 25:]; for reasons of completeness, in Appendix lAl we present its construction. In our simulations, the 
box is in the interval [1.57T, 2.57r], i.e., ipg.s.( x ii X2, x^) is zero whenever any x$ is outside of this interval. At t = the 
box potential is suddenly turned off, and the wave packet is simultaneously (and suddenly) imparted some momentum 
of magnitude K > towards the wall: ip(xi, x%, X3, t = 0) = ip g . s . (x±, X2, X3) exp[— iK(xi + X2 + £3)]; apparently, K 
denotes the imparted momentum per particle. From such an initial state, we are able to find projection coefficients 
&(fci, k.2, ks) defined in Eq. I|8[). that is, we can find the corresponding function G(fci, k2, ^3) which is needed to calculate 
the Fourier transform (|22p . The Fourier integral in (|22[) is in this particular example 3-dimensional, and it is calculated 
numerically by using the fast Fourier transform algorithm in MATLAB. This provides us with the time-dependent 
wave function ip(xi, X2, £3, t), which we use to study dynamics of observables such as the single-particle (SP) density 
p(x,t) or the momentum distribution n(k,t). 

First, let us explore the effect of the interactions on the reflections of a few-body Lieb-Liniger wave packet. In 
Figures [T] and [2] we plot the time-evolution of single-particle densities and distributions of the momenta, respectively. 
The plots are made at three different times, t = 0, 1, and 2, and for three values of the coupling parameter, c = 0.25, 
3, and 10. The magnitude of the imparted momentum per particle is K = 1. Note that the wave packets broaden 
in time due to the repulsive interactions between the particles, and also due to the wave dispersion effects; the wave 



7 





FIG. 3: (color online) Density evolution of a Lieb-Liniger wave packet comprised of N = 3 bosons, which is given some 
momentum kick K (per particle) towards the wall. Insets correspond to the times (a) t = 0, (b) t = 1, (c) t = 2, and (d) t = 3. 
The interaction strength is c = 1. Blue dashed-lines are for K — 0, red dot-dashed lines are for K = 3, and black solid lines 
are for K = 5. 



packets for larger values of c spread at a faster rate than the wave packets for smaller c. From Figs. Q] and [2] we 
observe that wave packets with a larger interaction parameter c get reflected faster than the wave packets for smaller 
c; for wave packets with smaller repulsion between the particles (smaller c), the compression of the wave packet is 
stronger, and therefore reflection of the momenta occurs at a slower rate. We also observe that all wave packets 
exhibit interference fringes during the reflection process. However, we find the interference fringes to be deeper for 
smaller values of c, which follows from the fact that the wave packets for smaller c are more spatially coherent. This 
can be seen also from Fig. [2] which displays momentum distributions. The distribution n(k,t) for c = 0.25, at the 
largest time shown t = 2, has one strong well-defined peak (the one closest to zero), and several smaller peaks of the 
wave components with larger magnitude of the momentum [sec Fig. |2Ja)]. In contrast, for c = 10 this most dominant 
peak close to k = is much smaller [see Fig. |U[c)]. 

Next we explore dependence of the time-evolution on the imparted momentum. To this end we fix the interaction 
strength at c = 1, and observe the time-evolution for three different initial conditions (see Figs. Q5]aiicl|4]) : (i) expansion 
in the presence of the wall occurs when K = 0, (ii) reflection at an intermediate value K — 3, and (iii) for large value 
of the imparted momentum K = 5. The wave packets for K = 3 and 5 have the property that basically all of the 
initial momentum distribution is directed towards the wall, i.e., the distributions at t = is on the negative fc-axis. 
In contrast, exactly half of the initial momentum distribution of the wave packet for K = is positive (negative). 
The basic distinction between these cases is that the wave packets with sufficiently large imparted momentum K get 
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FIG. 4: (color online) The same as in Fig. [3] but for the momentum distribution 



simply reflected from the wall and at larger times the interference fringes are almost negligible. For example, the wave 
packet with K = 5 is practically completely reflected from the wall at t = 2, see solid black lines in Figs. [3Jc) and0Jc); 
the momentum distribution is on the positive fc-axis and the interference fringes are essentially absent. In contrast, 
for K = half of the momentum distribution is already positive (corresponding to motion away from the wall), and 
this part interferes with the reflected component at all times of the evolution. Note that the wave packet with K = 
is still in the process of reflection from the wall at t = 2 because a large fraction of its momentum distribution is 
still on the negative fc-axis, see dashed blue line in Fig. Htc); the interference fringes are the largest in this case, see 
dashed blue line in Fig. EIc). 

Exact solutions can serve as a benchmark to check the range of validity of other methods which may be used to 
analyze noncquilibrium dynamics of interacting systems. We have compared the solutions obtained with the Fourier 
transform method presented here with the so-called hydrodynamic formalism (To| . which describes the Lieb-Liniger 
system via the nonlinear Schrodinger equation with variable nonlinearity |4lj |. In Figs. O (a)-(d), we show density 
profiles for two different couplings (c = 0.25 and c = 3) at two different times (t = 1 and t = 2). We find that the 
single-particle density (and momentum distribution), calculated within this method, are in good agreement with our 
simulations for small values of the coupling parameter c (up to c = 1); this upper limit for c also depends on the initial 
density of the ID Bose gas, as it is well known that the effective interaction strength parameter is c divided by the 
linear density [||. However, for larger values of c, the hydrodynamic formalism goes beyond its range of validity for 
the simulations presented here. For example, for the simulations at intermediate interaction strength c = 3 [sec Figs. 
[5] (c) and (d)], the hydrodynamic formalism predicts deeper interference fringes than those obtained via the Fourier 
transform method; this is attributed to the fact that the hydrodynamic formalism overestimates the spatial coherence 
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FIG. 5: (color online) Comparison of the density evolution in the exact calculation, with the hydrodynamic approximation 
(HDA) [(a)-(d)], and the Fermi-Bose mapping [(e) and (f)] valid in the Tonks-Girardeau (TG) regime. The interaction strengths 
c and times t in the insets are: (a) c = 0.25, t = 1; (b) c = 0.25, t = 2; (c) c = 3, t = 1; (d) c = 3, t = 2; (c) c = 10, t = 0; (d) 
c = 10, t = 1. The initially imparted momentum is K = 1 for all figures. 
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of the wave packet jlOj, |41|. For sufficiently large c, the system is in the Tonks-Girardcau regime, and one can employ 
the Fermi-Bose mapping fj\, l35j to study the dynamics. In Fig. [5] (e) and (f) we compare our calculation with that 
obtained via Fermi-Bose mapping (c = oo, 0,11a]) for a large value of the interaction strength c = 10; we observe that 
qualitative features of the Tonks-Girardeau regime such as the TV peaks in the initial single-particle density coincide 
in the two calculations, however, even larger c is needed to obtain better quantitative agreement. 

A. Normalization of eigenstates 

In order to numerically check our conjecture that the Lieb-Linigcr hard-wall eigenstates defined in ([5]) arc properly 
normalized, we have compared the initial state obtained via ip g , s .{x^i x 2i x s) ex P[ — iK(x\ + X2 + £3)], and the wave 
function obtained via Eq. (|12j) by employing the function G(fci,fc2,fc3), which is calculated from the projection 
coefficients b{k\, fe, k 3 ) as in Eq. (fl^|) . We found that the relative agreement between the two wave functions is on 
the order of 1% or better, which is on the order of the numerical accuracy for the size of our numerical grid, which 
is limited by computer memory. We have performed this comparison for various initial conditions (different K and 
c values). Unfortunately, a rigorous proof of normalization of Lieb-Liniger hard- wall eigenstates is to the best of our 
knowledge still lacking. 

B. Connection to physical units 

If we consider a system of 87 Rb atoms, then the ratio Xq/Tq « 3.65 x 10 _10 m 2 /s is fixed. By choosing for example 
X n w 1.35 /im for the spatial scale, the temporal scale is set to T = 5ms. The 3D scattering length is a = 5.3nm. The 
interaction parameter c can be varied by changing the width of transversal confinement for example, the values 
of c = 0.25 up to c — 10, can be obtained by varying l± from 242nm down to l± w 41nm, respectively. Of course, 
for a different choice of temporal and spatial scales, transversal confinements lj_ would have different values. We have 
also verified that for the choice of scales in our example, the longitudinal energy Eq is less then the transverse energy 
spacing huj±, a condition needed for freezing the radial degrees of freedom. 

V. CONCLUSION 

We have studied reflections of a Lieb-Liniger wave packet from the hard-wall potential. By employing the symmetry 
of the many-body eigenstates with respect to the change of the sign and permutation of their quantum numbers (i.e., 
quasimomenta) , that is, Equations ^ and (|10[) . we have demonstrated that time-evolution of this interacting many- 
body wave packet can be represented in terms of an TV-dimensional Fourier transform, where N is the number of 
particles in the wave packet. This result simplifies our understanding of the time-evolution in this many-body problem 
and enables straightforward calculation of the time-asymptotic properties of the system. 

We have utilized the formalism to numerically study dynamics of single-particle density and momentum distribution 
of a few-body wave packet reflecting from the wall (the wave packet is initially close to the wall). Reflection dynamics 
and interference phenomena depend on the strength of the interaction between the particles c and the imparted 
momentum K towards the wall. The wave packets for smaller c get reflected at a slower rate, because they get 
compressed more strongly as the wave packet hits the wall. Moreover, the interference fringes are deeper (larger 
visibility) for smaller values of c. If K is sufficiently large such that the initial momentum distribution is on the 
negative fc-axis, the wave packet gets reflected and the interference fringes become small as soon as most of the 
momenta become positive. On the other hand, for K = 0, the interference effects are fairly large. 
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APPENDIX A: THE GROUND STATE OF A LIEB-LINIGER GAS IN AN INFINITELY DEEP BOX [25J 

In section HVl we study Lieb-Liniger dynamics in the presence of the hard wall potential, with an example of three 
Lieb-Liniger bosons which are at t < confined in the ground state of an infinitely deep box of length L = ir. The 
ground state in fundamental permutation sector Ri has been constructed by Gaudin in Ref. (25j via a superposition 
of 2 N free space eigenstates. For the box in the interval [1.57T, 2.5w], the ground state (up to a normalization constant) 
reads 

■ • • cx £ * - * II f 1 - zrr) E II (i + ^ 9P ^- 1 - 5n) - ( A1 ) 

Here, summations are taken over 2 N elements of set {e}, and N\ permutations P. The quasimomenta qj = £j\qj\, for 
j = 1 , . . . , N, are determined by set of transcendental equations 



h\L = 77 + J2[t^- 1 — — +tan- 1 . (A2) 

1*1 - ml ml + ml ' 



Eqs. (| A2[) arc solved numerically. For the initial state corresponding to three particles (N = 3), where ip(xi, X2, %3, t = 
0) = ijj g , s ,(xi, X2, X3) exp[—iK(xi + X2 + £3)], it is straightforward to obtain the projection coefficients 6({fc}) by 
employing Eq. ||5J) and the orthonormality of eigenstates 4>{k}- 
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